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<Q Abstract 

in 



We consider the classical time evolution of a real scalar field in 2 dimensional 
Minkowski space with a A</> 4 interaction. We compute the spatial and temporal two- 
point correlation functions and extract the renormalized mass of the interacting theory. 



We find our results are consistent with the one- and two-loop quantum computation. 



We also perform Monte Carlo simulations of the quantum theory and conclude that 
the classical scheme is able to produce more accurate results with a fraction of the 
\Q ■ CPU time. 

o 
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1 Introduction 

In quantum field theory, the vacuum is filled with zero-point energy to which each mode 
contributes ftw/2. This picture raises the question, what dynamics would emerge in an 
interacting classical field theory if each mode is initialized with energy ftw/2. The dynamics 
of modes scattering and re-scattering off one another resembles the corresponding picture 
in quantum field theory. If we assume the phase of each mode is randomly distributed, 
then Lorentz invariant correlation functions arise upon averaging over phases. For example 
if we consider A0 4 theory, the value of <f? after averaging over phases equals the vacuum 
expectation value (0\(fi 2 (x)\0) of quantum field theory. This vacuum expectation value in 
combination with the A0 4 term produces a mass correction that is exactly the one-loop 
correction in quantum field theory. 

Classical zero-point fluctuations of free massless scalar and vector fields, the Lorentz 
invariance thereof, and the effects of phase averaging, have been considered before, most 
notably by Boyer pQ. By studying correlation functions of random electromagnetic fields, 
and the behavior of charged systems in their presence, Boyer has emphasized interesting 
similarities to quantum mechanics. 
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In an accompanying paper, two of the present authors have considered the dynamics of 
background fluctuations with an ftjjj/2 energy spectrum in interacting classical field theory. 
They showed that a lattice formulation allows for a nonperturbative analysis of the classical 
evolution equations. They found critical, strong coupling behavior in A0 4 theory that is 
strikingly similar to what is known about the quantum field theory. They also developed a 
perturbative expansion of the classical theory in the continuum and found a loop expansion 
very similar to, but apparently not identical to, that of quantum field theory. All of this 
motivates the present work, where we attempt a more quantitative comparison of the classical 
and quantum theories, where both are defined on the lattice. 

In the classical theory it is easy to implement the time evolution of the system, and we 
expect the configuration at a later time to incorporate the effects of interactions. From the 
2-point correlator we can obtain the physical mass, to be compared with the physical mass 
obtained by conventional methods from the lattice quantum theory. Our interest here is not 
in the continuum limit of these lattice theories, or the corresponding critical behavior. We 
wish to compare physical masses that are safely above the inverse size of the systems, and 
safely below the inverse lattice spacing. Our focus is on the loop corrections to these masses. 
This means that effects at two-loops or higher will need to be extracted, since we already 
know that the one-loop effects of the two theories agree. 

There are extensive studies of classical field theories and their simulation in the literature, 
typically to approximate quantum field theories in thermal backgrounds or other backgrounds 
with high occupation number El II]- In some of these studies one component of the 
initialized classical field does model the zero point spectrum of the quantum vacuum [Hj. 
But even then the focus has been on the dynamics of the thermal component, and not on 
how the interactions are manifested on the vacuum part. One of the common objects of 
study is the plasmon (or Landau) damping rate, related to the imaginary part of the on- 
shell self-energy j2|. Instead we compare to zero temperature field theory and study how 
closely the classical theory simulates quantum mass renormalizations, and thus our interest 
is in the real part of the self-energy. 



2 Time Evolution 

We study a classical A0 4 theory in 1+1 Minkowski space on a lattice. The action is 

s = E g _ WU + 1) - ««■ J))' _ ^d,j) _ (i) 

i j=0 

where we employ the units h — c — a — 1, and a is the lattice spacing along the x direction. 
The integers i and j label the sites in the t and x directions, and j runs from to N — 1. 
A periodic boundary condition is imposed along the x direction, i.e. (f>(i,N) = <p(i,0). 
The equations of motion are obtained from the action and we adopt a leapfrog method to 
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numerically integrate them. 

<f>{i + l,j) = 0(z,i)+a t 0(i + l/2,i), (2) 
<j>{i + 3/2, j) = < j>(i + l/2J) + a t ^>{i + l,j + l)-2 ( j>{i + l,j) + ( f>{i+l,j-l) 

-m^(i+l,j)-A0 3 (i + l,j)) (3) 

The leapflog scheme can be shown to be reversible and has second order accuracy. We 
obtain accurate classical evolution with a t = 0.1, where a t is the lattice spacing along the 
time direction. (Decreasing at didn't significantly change our results.) 

In quantum field theory, each mode contributes Huj/2 to the zero-point energy and thus 
the initial configuration of interest is 

^ 1 2nk 
MhJ) = 2^ rjr— cos(uj k a t i + — -j + fl fc ), (4) 

k=-N/2 V iVU;fc 



u k = J 4 sin 2 (vrA;/iV) + /x 2 , (5) 

where 9 k is a phase parameter, and u k is the lattice dispersion relation. We leave the mass 
/i as a parameter since the physical mass will be different from the bare mass m . One can 
easily check that the energy of each mode equals hw k /2 (with A = and fi = mo). 

The phases 9 k are independent parameters and we assume each 9 k is randomly distributed 
over [0, 27r). Then, for example, we obtain the following average values 

/d6 f d9' 1 
— - / — ^ cos fc cos 0£ = -6 kk >. (6) 
ZTT J Z7T L 

We shall use (•) to denote this average over the phases. This averaging gives an inter- 
esting indication of the emergence of quantum effects if we compute the average values of 

N ^ 1 2nk 
(0o(i,j)0o(i',j')> = 22 ^7 cos (^ a *( i + -/)) (7) 

k =-N/2 k 

= Re{0\T<f>(i,jW,j')\0), (8) 

where the r.h.s of the second line is the real part of the 2-point function in lattice quantum 
field theory, i.e. the real part of the Feynman propagator. This implies that the mass shift 
induced from the interaction term of the classical action is exactly the same as the one-loop 
mass correction in quantum field theory. 

We shall use (@J as the starting point for the investigation of the interacting theory. The 
key point is that the classical evolution of classical fields can be easily carried out in the 
presence of interactions and the configuration at a later time incorporates the interacting 
dynamics. We can use (HJ) to specify the initial conditions, 0(0, j) = (f> o (0,j) and 0(1/2, j) = 
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0o(l/2, j), but evolve forward in time with the full A 7^ equations of motion (J2J) and ©• 
We then compute correlators 



(9) 



where i and i' are sufficiently far away from the initial time. We obtain the average values 
by repeatedly computing </>(i,j)^>(i',j') starting from different sets of random phases. To 
reduce statistical noise we make use of translational invariance to also average over space. 

Various correlation functions can be calculated with respect to the final time slice of the 
simulation, at tf = a t ij. The time correlator G t (i), the space correlator G x (j), and the 
zero-mode time correlator Gq{i) are defined as follows. 



Any of these correlators can be used to extract mass, by fitting to the corresponding free 
field quantities. 

For an example we use N = 256, ml = 51/iV 2 , A = 20/N 2 , and t f = 3JV/4. We 
generate 5000 trajectories with an initial mass parameter /i = m , and obtain the correlation 
functions at tf. A least squares fit to free correlators with an adjustable mass produces a 
more appropriate input mass [i. We produce another 5000 trajectories with this fi as input, 
and then obtain a new /1. This process is repeated until the input mass agrees with the 
extracted mass. We find that it converges fast, typically in three interations. In this case 
the iteration converges to the value m? = 100. 0/N 2 . All three correlators have excellent 
fits to the corresponding free correlators with this mass. This is shown in Figs. (PUJ). For 
comparison we have also displayed the significantly different free correlators with /i = mo, 
to clearly show the effect of the mass renormalization. 

We shall find that almost all of this mass shift, to within a percent or two, is due to the 
one-loop renormalization. We will have to go beyond one-loop in the following. However 
this example has already demonstrated some important features of the classical simulation. 

• The system tends to evolve towards the renormalized physical mass from whatever 
mass is used to specify the initial condition. 

• The masses from the various correlators agree. This is an indication that Lorentz 
symmetry breaking due to lattice effects is not significant. 

• The physical mass can be determined to high accuracy. 
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Figure 1: The zero-mode correlation function versus position in lattice units. The 

zero-mode correlator from the classical simulation for itiq = 51 /N 2 and A = 20 /N 2 . It is 
well fit by analytical free correlator with fi 2 = 100/iV 2 (the lines completely overlap). Also 
shown is the free correlator with fi 2 = tUq. 
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Figure 2: The time and space correlators versus position (or time) in lattice 
units. The time (wavy) and space (smooth) correlators are shown for the same parameters 
as Fig. (1). Again we show the free correlators for comparison. 
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3 A Precision Measurement of Mass 



The one-loop graph calculated on a lattice with discrete space and continuous time, after a 
Wick rotation, is given by 



U(m) 



M ^ 



dk' 



-E 

ON ^ 



N ^ J 2tt k' 2 + 4 sin 2 (irk /N) + m 2 
1 



(13) 



2N k V4sin 2 {7rk/N) + m 2 ' 
This one-loop mass renormalization effect can be used to define the one-loop gap equation 



ni 



Sap 



m 2 + 3An(m gap ). (14) 



By obtaining the self-consistent mass m gap from this equation we effectively sum up graphs 
with chains of bubbles, where each bubble represents II. From the gap mass we may define a 
dimensionless coupling g = A/m 2 ap , and the extent to which the gap equation can be trusted 
will be determined by this coupling. The gap equation provides a useful representation of 
the physics for quite a large region of the rn^-X plane. This is because the one-loop graph is 
the only (log) divergent graph in this 1 + 1 dimensional theory, and it thus dominates. 

For the previous example we have g = 0.2 and m 2 ap = 101. 7/N 2 , to be compared with 
the mass extracted from the simulation of m 2 = 100/iV 2 . 1 Given that there are two- loop 
corrections to the gap mass that are not included in the one-loop gap equation, our goal 
then is to decide whether the difference between our extracted mass and the gap mass is 
consistent with the size of the two-loop corrections from quantum field theory. We will first 
need to discuss the errors and uncertainties inherent in the classical simulation, and work 
to reduce them, in order to obtain the required accuracy. Then we need to determine what 
the quantum corrections actually are, and for this we will consider both lattice Monte Carlo 
methods and lattice perturbation theory. 

At least at weak coupling we have seen that the evolving classical field configuration 
quickly settles down to a stable configuration that incorporates the effects of interactions. 
The stability of an interacting system with a nonthermal energy spectrum (an frw spectrum) 
is itself nontrivial, but at least in the continuum the frw spectrum is at a Lorentz invariant 
symmetry point, and thus stability may be a consequence. But on a lattice there is Lorentz 
symmetry breaking induced by the lattice itself, and this is expected to lead to thermaliza- 
tion. Thermalization of classical fields on a lattice is a well studied phenomena jl]. And 
consistent with these studies, our results indicate that the thermalization time scale is very 
long compared to the time scales of interest to us; on these shorter time scales, and for weak 
to moderate coupling, the interacting quantum-like configurations are very nearly free of 



1 Since m gap is much closer to m than is mo, in practice we use m gap as the initial trial mass [i to start 
the simulation. 
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thermalization effects. Presumably the lattice Lorentz symmetry breaking does not strongly 
affect the physics below the lattice cutoff. 

Given the observed stability at weak coupling we can ask how far we should evolve 
forward in time to obtain the correlators. We are motivated to average over a fairly long 
range of time for the following reason. When we extract the mass for different choices of tf 
we find that the extracted mass fluctuates slightly as a function of tf. For the example above 
it fluctuates within a ±0.5% band. This indicates that the configuration is not completely 
stationary, even if it is stable. But since the mass fluctuates around a stable central value, 
we may average the masses obtained over a range of tf and thus obtain an accurate mass. 

There is more of a problem for large couplings. Here we find that the central value can 
slowly drift in one direction as tf increases to large values, in which case we are unable to 
obtain a reliable mass. This could be taken as a sign of mild thermalization effects. 2 It is 
also sensitive to the size of the mass; for masses, m 2 N 2 > 80, the drift for increasing tf is 
towards higher mass, while for m 2 N 2 < 80 it is towards lower mass. 3 For m 2 N 2 up to a few 
hundred the drift doesn't show up until g > 0.6 

We have also seen that we can extract mass from three different correlators. Of the three 
we find that the space correlator is the least reliable. It shows the most sensitivity to numeri- 
cal errors and tf dependence. This is related to the fact that it is an exponentially decreasing 
function, and thus a fit is sensitive to only a small spatial range inversely proportional to the 
mass. The zero-mode correlator Gq{i) and the time correlator Gt{i) tend to agree well with 
each other. But the pure sinusoidal oscillations exhibited by Go(i) for small couplings tend 
to become damped oscillations at large coupling, due to the mild thermalization effects. 

We thus concentrate on the time correlator G t (i). Our prescription is to average over 
twenty equally spaced values of tf up to a final time tf = ION (a temporal extent 10 times 
the spatial extent) to obtain a mass. We then average between 10iV and 20iV to obtain a 
second mass. If the second mass does not deviate significantly from the first, then we accept 
the first mass as a reliable extracted mass. We also note that the time correlator contains 
useful information in its deviation in functional form from the best fit free correlator, this 
effect becoming visible at larger coupling (it is not visible in the weak coupling correlators 
in Figs. (PEJ))- These deviations in form are expected since the mass corrections at higher 
loops are momentum dependent. 

2 At couplings significantly larger than one these effects are more pronounced, and the amplitudes of the 
lowest modes, and in particular the zero- mode, can be seen to grow with tf. 

3 This is for intermediate size coupling; for couplings significantly larger than one the drift is towards 
larger mass. 
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4 Comparison to Quantum Methods 



We now turn to the extraction of the mass in a lattice quantum theory. We would like to 
use an action and a lattice as similar as possible to the classical theory. However, standard 
lattice methods for the quantum theory make it necessary to Wick rotate the time direction, 
and a periodic boundary condition is imposed along it as well, giving a periodic Euclidean 
space. Correspondingly we take Of — a — 1 rather than the 10at = a = 1 in the classical 
simulation. Then the action for lattice quantum field theory is 

s = j: m + hJt: gfrjff + WLi ± *Ml + aj^ hj) + a 04( . . } (15) 

The different lattice regularization causes an obvious difference in the value of the one- 
loop self-energy graph (to be compared with ([T3])). which now becomes 

nE ( m ) = w X) 4 sin 2 (vrfc /iV) + 4 sin 2 (ttA:' ~Jn) + m 2 ' ^ 
The gap equation for this theory (where the subscript E stands for the Euclidean lattice) is 

m lga P = m l + 3An E (m Egap ). (17) 

The resulting values m Egap and g E w iH thus differ from m gap and ^ from our previous gap 
equation. 

It is useful to graphically compare the results of these gap equations. Fig. (jSJ) and (JH) 
show the behavior of lines of constant g and g E in the m^-A plane. We see that large g occurs 
when m 2 , is sufficiently negative. The largest coupling line on this graph is g — 10.24, which 
is the value of the critical coupling as determined by previous Monte Carlo simulations of the 
quantum theory |S]. 4 In contrast to the gap mass, the physical mass is expected to vanish 
on this line. We see that the g E = 10.24 line has quite a different behavior at small A since it 
does not intersect the origin. The reason for the difference is that n E (m) — > 1/m 2 as m — » 
while n(m) — > 1/m. 

We further illustrate the differences in the two lattice regularizations in Fig. where we 
compare the lines of constant gap masses. We see that the differences basically grow linearly 
with A. We would like to account for this obvious source of difference when extracting the 
physical masses. We therefore concentrate on the deviation of the physical mass in each 
theory from the respective gap mass. The leading order deviation from the gap mass is 
a finite two- loop effect, and we expect that the lattice dependence of this leading order 
difference will be small. This gives us reason to focus our study of the physical masses at 
weaker coupling. 



4 The critical coupling obtained from density matrix renormalization group methods 6 is 9.98. 
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Figure 3: Lines of constant g, with A and Wq in units of 1/N 2 (N = 256). The g of adjacent 
lines decreasing from left to right by a factor of 0.8. The largest coupling is the critical 
coupling of g = 10.24 and the smallest is g = 0.12. 
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Figure 5: Lines of constant tri 2 ap and ^ggap * n units of 1/N 2 . The mass values are given by 
the intersection of these lines with the x-axis. Also shown is a special constant mass line 
and the critical coupling lines of the two theories. 

Also displayed in Fig. © is a special line of constant gap mass, where the gap masses of 
the two theories agree. For iV = 256 this occurs at m^ ap = m| gap = 4.9338/iV 2 . And since 
the gap masses agree, g = along this line. Thus in the vicinity of this line the obvious 
differences between the two theories are minimized. In particular the two theories would 
agree on the location of the critical point along this line. 

In our lattice quantum simulations we start with a hot start, and as our updating scheme 
we choose the heat-bath algorithm followed by a Wolff step to reduce critical slowing down. 
'One update' of the lattice consists of visiting each site with the heat-bath algorithm and 
then applying a Wolff single-cluster step [7|. We monitor the relative size of the generated 
Wolff clusters, and we ensure that the system has thermalized by monitoring the action. 
The autocorrelation time r is estimated by measuring the autocorrelation function p(t) of 
various operators, such as 0(0, 0)0(0, N/2) and the action, and computing the integrated 
autocorrelation time according to T int — \ + J2t=iP(^)- Pi®) = 1 an< ^ M is such that for 
t > M, p(t) is essentially noise. 

The physical mass is obtained by fitting the correlation function 



to a cosh or a decaying exponential. Here (■) is the thermal expectation value over different 
configurations. The systematic errors are the most difficult to control. To monitor these we 
compare results at neighboring values for parameters where similar results are expected. We 




(18) 
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also repeat the same simulations with different random number generator seeds. 

As a check on the quantum simulation we find that it does give physical masses compatible 
with the gap equation, to the extent that the gap equation can be trusted. The mass 
extracted at strong coupling is also consistent with the expected critical behavior. The 
decrease of the physical mass towards zero as is increased serves as one method for the 
determination of the location of the critical line. For another, the order parameter 




i,j=0 

and its histogram can be studied as the size of the lattice is changed (see [5 ). The same 
data can be used to extract a critical exponent v defined by 

mp h y S oc |A - A c | y . (19) 

In all these studies we obtain results from our quantum simulation consistent with previous 
studies. 

But the quantum simulation proved to be rather deficient at producing the type of pre- 
cision that is needed in our study. A meaningful comparison to the classical theory requires 
that we accurately determine the difference between the physical mass and the gap mass. 
For a 256 x 256 lattice, and for a run time of the order of a few days to extract a mass 
at a single (mo, A) point, a large percentage error still arose in the mass difference. 5 This 
error seems to be related to problems with correlated data and/or autocorrelation, and it 
is exacerbated for smaller physical mass. This ruled out any meaningful comparison of the 
theories along the special line m^ ap = m| gap = 4.9338/iV 2 that we identified above. But as 
we shall see shortly, the quantum simulation results will still serve a useful purpose. 

5 The Two-Loop Effect 

We turn to improving the gap equation itself. Although the gap equation is nonperturbative 
in the sense that it resums an infinite set of graphs, it omits some effects of second order 
and higher in the coupling. At leading perturbative order the difference between the physi- 
cal mass and the gap mass is determined by a single two-loop diagram, the "sunset" diagram: 




Through the use of a multi-processor cluster we have looked at roughly 50 points, using 50000 configu- 
rations for each. 
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In the lattice version of the quantum theory, on the square Euclidean lattice, this diagram 
is given by an explicit multiple summation, 

N/2-1 

n 2 E{m,k) = — G(Pi + k,P2)G(q 1 ,q 2 )G(p 1 + q 1 ,p2 + q 2 ), (20) 

{pi,P2,qi,q2} = -N/2 



G(pi,p 2 ) = — - t j^r—r-^T—^TT—,- (21) 



1 

4 sin 2 (7r^i/iV) + 4 sin 2 (7rp 2 /A0 + m 2 ' 

Here k represents an external momentum. This calculation gives a correction to the 2-point 
function at what corresponds to space-like external momentum, since this is a Euclidean 
quantum theory, and we find that n 2 E(w, k) monotonically increases as k —>■ 0. The classical 
theory however has Lorentzian signature, and so we actually need II2 at the time-like on- 
shell momentum. We will use n^E^, 0) as an estimate for this. 6 Moreover our calculation 
is on the square Euclidean lattice, rather than the lattice of the classical simulation, but this 
additional error should be small. 

For a g 2 improved gap equation to be useful we need to have some idea of the range of g 
over which it is accurate, before even higher order effects become important. This is where 
the quantum lattice results are useful. Although there are large errors on the extracted 
mass correction at any single point in parameter space, when taken collectively the quantum 
simulation data supports n 2 E(^, 0) as an estimate of the difference between the physical 
and gap masses. This appears to be a good estimate for couplings as large as g ~ 0.6. 
As a further check we have directly estimated the single three-loop graph using the same 
technique as above (for smaller lattices) and found that it is < g/2 times the size of the 
two-loop graph. 

We therefore consider a two-loop gap equation of the following form, 

m ga P = m l + 1 1 - £("v P )]3n(m gap )A - 6n 2E (m gap , °)^ 2 - ( 22 ) 

We have inserted a correction factor in the A term to account for possible lattice artifacts 
or finite size effects inherent in our classical simulation. 7 It is important to account for this 
because a small correction in the A term can compete in size with the A 2 term. Similarly 
a correction to the A 2 term of similar size can safely be ignored. This gap equation can be 
considered as a quantum model to be tested against the classical simulation data. 

Our strategy then is to use ()22)1 as a one parameter model for the extracted physical mass 
m p h ys from the classical simulation. We run the simulation for a range of values of tuq and 
A that produce values for m p h ys very close to a fixed m gap , and thereby determine e(m gap ) 
through a best fit. We choose the couplings g = (0.1,0.2,0.3,0.4,0.5). We use 5000 trajec- 
tories for each determination of a mass. For the four values ^ gap = (60, 80, 100, 140)/iV 2 we 
find£(m gap ) « (0.032,0.029,0.029,0.029). 

6 By naive extrapolation from space-like momenta it is an underestimate, but probably by less than 10%. 
7 In particular there is the effect responsible for the slightly fluctuating nature of the mass as a function 
of tf that we have noted earlier. 
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To test the quantum model for m p h ys we write it in terms of the dimensionless coupling 
g = A/m 2 ap and then isolate the g 2 dependence to get 

g2 = m ph ys - { m l + [1 - g(^ g a P )]3n(m gap )m^ p ^} ^ 
-6n 2E (m ga p,0)m^ ap 

We then determine the r.h.s from the classical simulation for a fixed m gap and for a range 
of m and g that give m p h ys ~ rn gap . This is displayed in the plots in Fig. © where we 
see the residual quadratic dependence on g, and compare to the superimposed g 2 line of the 
quantum model. (The points would lie along a straight line if there was no g 2 effect.) 

The spread of the points for 10 different random number seeds gives an indication of the 
errors. Of the four choices of m gap , the errors in the m 2 &p = 80 /N 2 (60 /N 2 ) plot are the 
smallest (largest), and for larger masses the errors again grow larger. In section 3 we have 
already noted that m 2 ap = 80 /N 2 suffers least from the mass drift phenomenon at strong 
coupling, and thus the pattern of errors observed here are correlated with that phenomenon. 
In fact the drift is found to occur for g > 0.4 in the m 2 ap = 60/iV 2 case, 8 and we can see 
how the sign of the drift tends to skew the g 2 effect towards a larger value. 

The other three plots show that the results of the classical simulation are very well 
described by the quantum model, where the coefficient of the A 2 term is given by the sunset 
diagram of quantum field theory. This is our main result. We can also consider a two 
parameter fit to the data in these three plots, where we allow the coefficients of both the A 
and A 2 terms to vary. This yields a A 2 coefficient ~ 8% larger than the quantum prediction. 
For an independent set of data for masses m gap iV 2 from 60 to 200, where we only considered 
values of g up to 0.3, we obtained values for the A 2 coefficient that deviated between —2% 
and 13% from the quantum value. All these deviations are similar to the theoretical error 
in our estimate of the two-loop graph on mass shell. 

It is instructive to consider again lines of constant m gap on the rn^-X plane, where we 
now use the corrected two- loop gap equation (J22J- In Fig. (jJJ) we compare the result to the 
original one-loop gap equation, by showing the effect of the e correction with no A 2 term, 
and then the effect with the A 2 term added as well. This figure highlights the fact that the 
simulation has to be highly accurate to see the small A 2 effect. We have attained this level of 
precision in the classical simulation, and we have come nowhere close to it using conventional 
Monte Carlo methods in the quantum theory. 



6 Summary 

We have discussed the extraction of quantum loop effects from a classical simulation. We 
considered the time evolution of a scalar field configuration realizing the zero-point energy 

8 According to our previous prescription in section 3 the g = 0.4 and 0.5 points should not be considered 
reliable. 



13 



0.3- 
0.25- 


60 


e 


0.25- 




80 


0,2- 
0.15- 

0.1- 
0.05- 






0.2- 
0.15 - 

0.1- 
0.05 - 




°v 





0.1 oi 


0.3 0.4 0.5 





i 

, — -. : a — _ 

0.1 


0.2 0.3 0.4 0.5 


0.25- 






0.25 - 






0.2- 
0.15- 

0.1- 
0.05- 


100 


A 


0.2- 
0.15 - 

0.1- 
0.05 - 



1 — — , — - — r"T~~r 


140 : / 

8 

B 





0!l 02 03 04 0^5 


0)1 0.2 0.3 0.4 0.5 



Figure 6: Plots of the r.h.s of Equation (1231) versus g showing a comparison 
between the classical simulation and the quantum model. The graphs are labelled 
by m 2 N 2 and the solid line is the g 2 quantum field theory prediction. Each plot displays the 
results of 10 different random number seeds. 
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Figure 7: Lines of constant mass, for g up to 0.5. The bottom line in each case comes from 
the original gap equation, the next includes the e correction, and the top includes the A 2 
correction as well. 

spectrum of quantum field theory and including an averaging over phases. The time evolu- 
tion incorporates the effects of interactions and the trajectory quickly approaches a stable 
configuration on the time scale of interest, at weak to moderate couplings. The physical 
masses obtained from either the time or space correlation functions are in agreement, as 
would be expected by Lorentz invariance. Our method is clearly distinct from quantum 
simulations and uses only a tiny fraction of CPU time. This fraction is roughly 10~ 4 . 9 

In this way we have been able to perform a precision measurement of the effects of two- 
loop corrections within a classical simulation, with results that are consistent with quantum 
field theory. In jHj two of us studied the classical simulation at strong coupling and looked 
at two signatures of critical behavior, the overlap of the space and time correlators signaling 
masslessness and the development of a vacuum expectation value for the field. The findings 
were consistent with the location of the critical line in quantum field theory. 

On the other hand in jHj it was also indicated that the loop expansion of the classical 
theory was not identical to the quantum theory; and it was suggested that the sunset dia- 
gram, the particular two- loop diagram of interest here, may have missing contributions in the 
classical theory. Due to uncertainties our results here may still be compatible with missing 
contributions, but they imply that the missing contributions in this case are accidentally 
small or vanishing. More generally, any perturbative discrepancy is not easy to reconcile 
with the success of the classical simulation at both strong and weak coupling. However this 

9 We compare times to produce comparable accuracy, 500 trajectories in the classical simulation versus 
50000 configurations in the quantum simulation. 
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puzzle is resolved, we still find the concept of a classical simulation recovering most or all of 
a two-loop effect to be rather surprising and unexpected. 

Note added: For further work on resolving this puzzle see JH]- 
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